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Nomenclature 


AR 

* 

c,c 

C a , C L 

C D 

C Di 

C D0, C f 

Cp 

Cq 

Cd s 

Cd w 

D 

e 

F 

h 

k 

m 

M 

P 

Ps 

Pm 

PSFC 

q 

R 

Re, R 
Rf 

s 

Si 

SFC 

T 

TSFC 

u 

U 

V 


wing aspect ratio 
chord length 
coefficient of lift 
coefficient of total drag 
coefficient of induced drag 
coefficient of parasite drag 
coefficient of pressure 
volume coefficient of suction 
coefficient of suction drag 

coefficient of wing drag 

airframe drag 
lift efficiency 

dimensionless frequency, thrust 

heating value of fuel 

1/pARe 

mass flow rate 

Mach number 

power, pressure 

total ideal suction power 

power required to maintain flow 
power specific fuel consumption 
dynamic pressure 
range 

Reynolds number 

recovery factor 

wingplanform area 

equivalent suction area 

specific fuel consumption * 

temperature 

thrust specific fuel consumption 

x- component of the velocity in the boundary layer 

velocity in the freestream 

velocity 


ill 



W weight 

We airplane fuel empty weight - 

Wq airplane take-off gross weight 

Greek Symbols 

b flap deflection angle 

s spatial amplification rate of disturbances 

r density 

t thrust 

w, u suction velocity normal to wing surface area 

y angle of wave number vector 

h efficiency 

D incremental change 


subscripts 

A 

bl 

C 

d 

e 

eff. 

f 

i 

LFC 

m 

0 ,- 

P 

s 

ps 

t 

tr 


airplane A 

boundary layer thruster 
airplane C 
total drag 
exit plane 

effective drag or drag coefficient 
fuel 

inlet to propulsion system 
laminar flow control 
maximum 

overall or stagnation value, freestream 

propulsion 

suction 

pumping system 

thermal x 

transfer 
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1 . Introduction 

Reduction of drag has been a major goal of the aircraft industry as no other single 
quantity influences the operating costs of transport aircraft more than aerodynamic drag. It 
has been estimated that even modest reduction of frictional drag could reduce fuel costs by 
anywhere from 2 to 5%. Current research on boundary layer drag reduction (see Bushnell 
[1]) deals with various approaches to reduce turbulent skin friction drag as a means of 
improving aircraft performance. One of the techniques belonging to this category is 
laminar flow control in which extensive regions of laminar flow are maintained over aircraft 
surfaces by delaying transition to turbulence through the ingestion of boundary layer air. 
While problems of laminar flow control have been studied in some detail the prospect of 
improving the propulsion system of an aircraft by the use of ingested boundary layer air 
has received very little attention. 

An initial study for the purpose of reducing propulsion system requirements by 
utilizing the kinetic energy of boundary layer air was performed in the Mid-1970's at 
NASA Lewis [2]. This study which was based on ingesting the boundary layer air at a 
single location, did not yield any significant overall propulsion benefits and therefore the 
concept was not pursued further. However, since then it has been proposed that if the 
boundary layer air were ingested at various locations on the aircraft surface instead of just 
at one site an improvement in the propulsion system might be realized. The present report 
provides a review of laminar flow control by suction and focuses on the problems of 
reducing skin friction drag by maintaining extensive regions of laminar flow over the 
aircraft surfaces. In addition, it includes an evaluation of an aircraft propulsion system that 
is augmented by ingested boundary layer air. 

1.1 The LFC Concept 

The laminar flow control concept consists of making use of the available kinetic 
energy of ingested boundary layer air in augmenting the thrust provided by the main 
engines through an auxiliary device which we shall call a boundary layer thruster (BLT). 
The components of the boundary layer thruster are a suction compressor and an auxiliary 
nozzle. The suction compressor is powered by one of the turbines of the main air- 
breathing engines. The operation is much like that of an aft-fan operating by the input of 
power in a turbofan engine by a turbine. The transmission efficiency involved in the 
conversion of output power from the turbine to the suction compressor is considered to be 
100%. Boundary layer air is withdrawn from the upper and lower surfaces of each wing 
through suction slits into a collection surface as the schematics of Figure A.l indicate. The 
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air is then conveyed by suction ducts and compressed internally by the suction compressor 
before it passes through a converging-diverging nozzle ( in the case of supersonic speeds ) 
before being ejected aft of the fuselage to provide the additional thrust. This simplified 
thrusting system allows us to model the Y internal fluid mechanics in a simple way and to 
provide a measure of the performance of the thruster as a function of the exit velocity ratio. 
The suction slits are aligned in such a way to provide maximum flow per unit slot area and 
are sized and spaced spanwise so that they do not contribute to transition of the laminar 
boundary layer to a turbulent one ( see [9] ). The detailed construction of the collection 
surface has previously been considered by Boeing [9]. The relations involved in 
determining the system performance, the results obtained and the computer program used 
are provided in Appendix C. 
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2 . Review of the available literature on laminar flow control 

Under certain conditions undesirable dead regions are created in liquid and gas 
flows. These cause very prejudicial losses of energy. These losses can be avoided or 
reduced by drawing off small quantities of fluid from the surface into the interior of the 
body and thus preventing the development of turbulent regions. Boundary layer control by 
suction applied to wings insures an increase in maximum lift and permits using thick wing 
sections without excessive wing section or profile drag. One of the earliest experiments 
were carried by O. Shrenk [3] who investigated a large number of different arrangements 
of suction slits and their effect on maximum lift. It was shown by Shrenk that high lift 
coefficients were obtainable at certain suction volumes. In 1940 Shrenk [4] showed that 
boundary layer suction was more 'favorable' for airfoils with higher thickness ratios, in the 
sense that lower suction volumes are needed to achieve the desired lift coefficient (and 
hence a smaller expenditure of suction power) for an airfoil with a 20% thickness ratio in 
contrast to an airfoil with a 12% thickness ratio. Subsequent photographs taken during the 
flight test also show that an increase in the angle of deflection of the flapped airfoil is more 
favorable for laminarization of the upper surface of the airfoil and subsequent increase in 
lift and reduction of drag. However, the main thrust of the experiment was performed with 
the objective of studying the variation of lift with suction and this suggests to us that any 
study be it experimental or otherwise should be directed or aimed at maximizing certain 
aspects of aircraft performance. If range is the parameter of interest then the study should 
be carried out with the intent of maximizing the product of Bp /Cp for a conventional 

airplane. It should be noted that the studies [3] and [4] were carried out with suction 
employed on the upper surface of the airfoil in contrast to a study aimed primarily at the 
reduction of drag which would involve suction over both surfaces of the airfoil. Since in 
this case depending on the suction distribution there would be relatively no appreciable 
increase in maximum lift as acceleration of the boundary layer takes place on both surfaces 
with little relative (with respect to the upper surface) increase of the pressure distribution on 
the lower surface of the airfoil. One of the earliest flight tests that enhanced the prospects 
of laminar flow control by suction was the Miles experimental airplane. These tests 
indicated a reduction of about 22% in profile drag. 

Later in 1946 Smith and Roberts of the Douglas Aircraft Company published their 
findings [ 5] regarding the prospects of laminar flow control. In that paper they discuss the 
possibility of the separation of the main boundary layer flow for the potential flow over a 
flat plate obtained by the superposition of a rectilinear flow and a sink; in particular, that the 
increase in pressure downstream of the slot may result in an adverse pressure gradient and 
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subsequent separation. Thus, if the flow into a given slot is increased the closing 
streamline reaches further into the higher energy layers of the local flow field and a 
possibility of unstable flow appears. The authors mention that separation has been 
encountered by them in tests of boundary layer on a wing. It is surprising because it 
contradicts the usual belief that boundary layer removal always reduces separation. The 
increase in friction coefficient accompanying increased suction may be adequately 
demonstrated by Figure 1 of Appendix D, which however does not account for separation, 
as reported by Schlichting [6] from the results obtained by R. Iglish [7] for the continuous 
suction over a flat plate from a laminar boundary layer. It is interesting to note that with 
higher suction flow rates the friction coefficient becomes independent of viscosity and the 
drag obtained is that due to the sink effect of suction for a body immersed in a frictionless 
flow. 

Increased suction results in increased suction power and hence increase in total 
power expenditure in keeping the flow system in operation. This fact is the main cause that 
suction systems should be designed to operate at the minimum sufficient suction power 
needed to keep the flow stable. However, the flow pattern over a wing is more 
complicated than that over a flat plate and an optimum suction distribution is required to 
minimize disturbances and keep the flow pattern laminar. The above paper also mentions 
that the flying qualities of a jet aircraft with boundary layer removal are exceptionally good 
particularly during takeoff and landing. In relevance to this is the reduced angle of 
incidence for a flapped airfoil at a given flap deflection with boundary layer suction. It 
should be mentioned however that the title " The Jet Airplane Utilizing Boundary Layer Air 
for Propulsion " is somewhat misleading in the context of propulsive enhancements, for the 
paper is based in its entirety on aerodynamic performance and the reduction of drag and 
increase in maximum lift rather than the use of boundary layer air for propulsion. Stability 
analysis and suction flow rates are not given due to the lack of sufficient research pertaining 
to air requirements and stability of boundary layer flows at that time. It is interesting to 
note the results presented regarding the hypothetical study of the application of boundary 
layer control to actual aircraft reveal a considerable increase in range with the use of 
boundary layer suction. The performance comparison for three hypothetical airplanes are 
presented in Figures 2 and 3; airplane A with reciprocating engines, B with turbojet engines 
and C with turbojet engines with boundary layer inlets. It is found that the boundary layer 
control jet airplane excels the conventional in its payload carrying capability to about 2020 
miles, whereas a jet airplane with ramming intake can excel it up to only 1,310 miles. 
These ranges are carried at 430 m.p.h., 400 m.p.h., and 200 m.p.h., respectively, for the 
boundary layer intake jet, ramming intake jet, and conventional airplane. 
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Pfenninger in 1949 (see[8]) showed, with the use of suction slots, that with ingestion of 
small quantities of air (Cq= 0.0014 to 0.0018; where Cq= - [pco] s / [pu]oo) over both 

surfaces of an airfoil of 17% thickness ratio the boundary layer can be kept completely 
laminar with a reduction in the profile drag to about one-half its turbulent value at Re = 2.4 
x 10^ . Some of the problems faced in maintaining laminar boundary layers will be 
discussed briefly in the following introduction. 
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Flow Problems Over Swept wings 

The drag of an aircraft at cruise flight conditions is about 60% friction drag [ ... ]for 
present-day transport aircraft with turbulent boundary layers on their wetted surfaces. For 
underwater vehicles, the friction drag is about 90% of the total drag. In each case 
laminarizing the boundary layer offers substantial improvement in surface friction which in 
some cases may amount to 50% reduction of profile drag. 

Drag reductions of this magnitude are possible using extended natural laminar flow 
(NLF) or controlled laminar flow (LFC). The former principally applies to maintaining 
favorable pressure gradients along the wing surface to accelerate and stabilize the boundary 
layer by minimizing surface waves and discontinuities. The latter relies upon suction 
through slots in the wing or suction through a porous surface in order to prevent transition. 
The hybrid concept is another method which has evolved over the past years. Hybrid 
laminar flow control (HLFC) technique is a means of reducing airplane wing friction drag 
by combining suction laminar flow control near the leading edge ( forward of the front 
spar only ) with pressure distribution tailoring or natural laminar flow in the midchord 
section. This allows for maintaining laminar flow up to 75% wing chord. It appears that 
full chordwise suction inhibits the effectiveness of laminar flow control and lessens the 
benefits of drag reduction. On the other hand, sweep angles typical of modem commercial 
transport aircrafts are still somewhat higher than those for which substantial NLF has been 
demonstrated. This is how the hybrid concept originated. This method was successfully 
used by the Boeing Commercial Airplane Company (BCAC) in 1982 (see [9]). 
Considerable savings in wing profile drag with the deployement of HLFC can be seen from 
Figures 4 and 5 below. 

Basic to the use of suction in laminar flow control is an understanding of the 
mechanism of boundary layer transition and the types of instabilities that occur over a 
swept wing (see [10]), since stability of the flow dictates the maximum allowable suction 
flow rate that is to be ingested. However, as we shall see in 4.1 wherein the concept of a 
propulsive efficiency greater than one is introduced and the relative comparison between a 
boundary layer ingesting engine and an air-breathing engine is made, the concept of a 

’■"-v 

maximum value of suction flow is irrelevant. On a high speed swept-wing four basic types 
of boundary layer instabilities can occur: 1) viscous or " Tollmein- Schlichting " instability; 
2) inflectional or cross-flow instability; 3) dynamic or " Taylor-Goertler " instability; and 4) 
leading-edge attachment line contamination. 

Tollmien-Schlichting (T-S) instability depends upon the action of viscosity to 
transfer energy from the mean flow to the boundary disturbance. Amplification of T-S 
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disturbances is small in regions of favorable pressure gradients and large in regions of 
adverse pressure gradients. Moderate suction quantities may be employed to stabilize the 
T-S disturbances. 

The strong flow acceleration in the leading-edge region of swept wings induces a 
severe boundary layer cross flow and requires strong local suction. The effect of sweep 
and pressure gradient on the amplification ratio of disturbances is clearly seen from the 
results of Runyan and Steers [11]. 

Taylor-Goertler instability occurs primarily in the flow over concave surfaces. 
Most supercritical wing sections however do not contain any concave surfaces; as such, 
this type of instability is not a factor. 

Turbulence originating from a leading-edge roughness or a turbulent boundary layer 
starting at the wing-fuselage intersection may spread in both the spanwise and chordwise 
directions. This is referred to as leading-edge attachment line contamination and may be 
reduced substantially by reducing sweep angle and by removing the entire turbulent 
attachment- line boundary layer and re-establishing a laminar layer by means of suction. 

Particular attention is paid to the cross flow and T-S instabilities in connection with 
the results of Mack (see [12]) based on his investigation of the stability of the laminar 
boundary layer on two transonic wings of infinite span with distributed suction. Both 
wings have supercritical airfoil sections; one has a sweep angle of 23 deg with M^ = 0.82 

and c = 1.96m, the other has a 35 deg sweep with M^ = 0.891 and c = 2.0m . It is seen 
from Fig. 6 for the 35 deg wing that their is a considerable shift of the curve of maximum 
spatial amplification rate (a m ) of stationary (F = 0) or cross flow disturbances downwards. 
It also appears that with suction applied at the leading edge the peak of the curve decreases. 
The effect of the sweep angle on the cross flow is evident from the comparison of the 
results of the 35 deg wing with those of Fig. 7 for the 23 deg wing where the maximum 
value of <r m is less than 5 x 10'^. Compressibility effects tend to dampen the T-S 
disturbances (y = 0) in the mid-chord region. The design suction of the 23 deg wing is 60- 
70% of that of the 35 deg wing but it still controls the instability of these waves in this 
region as is shown by the large decrease in a m with suction. The pressure gradient on the 
upper surface of the 23 deg wing in the vicinity of the mid-chord region is nearly zero. 

So far the foregoing discussions have focused on the problems of laminar flow 
control and a review of some of the literature delineating the benefits in drag reduction 
offered by LFC. Nothing has been stated regarding the possible improvements in the 
propulsive system of an aircraft with LFC. The object of the following discussions will be 
to examine the boundary layer from a propulsion point of view and to show that the kinetic 
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energy of boundary layer air may be utilized in augmenting the propulsive efficiency of an 
aircraft by the use of a boundary layer thruster consisting of a suction compressor and a 
converging nozzle. The range equation for a non-conventional aircraft (one using 
boundary layer air for propulsion) will be* re-examined taking into effect the improvements 
in the propulsive efficiency and the L/D ratio . 
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4. Performance Prediction , 

4 . 1 The available energy in the wake of a boundary layer and its utilization in the 
increase in propulsive efficiency 

Ackeret [13] seems to have been the first to demonstrate the inherent advantage of 
the withdrawal of the boundary layer into the surface of a body which is propelled through 
fluid. He drew upon the fact that the available wake kinetic energy in a boundary layer is 
substantial. For simplicity consider a two-dimensional flat plate in a stream ( U,0,0) of 
incompressible flow as shown in Fig. 8. Let the velocity near the trailing edge be (u,v,0). 
The drag per unit span is given by, 

D = j pu(U-u)dy (1) 

and the power P m required to maintain the flow is. 

Pm = pU 3 cf <£)U-£)d(J) (2) 


If suction is now applied whereby the entire fluid within the boundary layer is withdrawn at 
the trailing edge into the surface and its total pressure is restored to the freestream value, 
which assumes no frictional losses in the ducting system and considers the pump efficiency 
to be 100%, then the total suction power required is the rate of change of kinetic energy 
between the inlet to the slots at the trailing edge and the exit to the duct. This suction power 
is given by 

P s =| pu 1/2 ( U 2 - u 2 ) dy (3) 


or 


P S = P U 3 C (JL) 


(1-iL) |l/2(l+JL)\d(I) 


m 


U 


(4) 


Comparison of eqs. (3) and (4) reveals that P s < P m . The term P s is also called the 
unavailable thermal energy ( see Ackeret [13] ). For ^ laminar boundary layer ( Blasius 
profile) this is about 78.7% of the product of plate friction drag times the freestream 
velocity and the available wake kinetic energy is about 21.3% of P m . It is thus seen that if 
the flow in the wake is ingested and ejected with a freestream velocity U the entire 
momentum loss in the boundary layer can be overcome and the suction power required to 
do so would be less than that needed to maintain the flow system. 
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The ratio of the thrust , obtained in acceleration of ingested flow from the 
boundary layer , to the frictional drag, may be shown using Fig. 1 by the following 
numerical examples. Assuming an asymptotic suction profile with C q = 3.75 x 10' 3 , C f 
= 0.9 x 10' 4 the value of the thrust to drag ratio is approximately 83.3% when the 
ingested flow is accelerated to freestream velocity. With a suction quantity of 14 x 10' 4 
the ratio is 99% when the suction flow is accelerated to a velocity of V e = 1.2 U^. 

Greater gains are obtainable when suction is applied steadily over the plate and 
laminar flow is maintained in the boundary layer instead of a turbulent flow. This point 
will be emphasized in subsequent discussions. It is also seen that the ratio P m / P s 
provides a measure of the propulsive efficiency which in this case is greater than one. The 
Froude efficiency of propulsion has always been defined as the ratio of the useful work 
done to the energy imparted to the fluid and the attainment of a value greater than one is not 
mere speculation especially in view of the principle of aircraft propulsion by wake 
regeneration as has been pointed out by Davidson [14] and its utilization in the Ogee 
scheme. 

4 . 2 The Range of a Subsonic Airplane With Actively Controlled Boundary Layer From 
a Propulsion Point of View 

The same bleed-off boundary layer air that is used to improve the lift-drag ratio of 
the airframe is considered as the inlet air of the auxiliary propulsor to increase its 
efficiency. The primary points made, that are in disagreement with some of the existing 
studies are: 

(1) that it is incorrect to charge the same friction drag of the airframe to the inlet 
momentum of the propulsor. 

(2) While it may give the correct range it does not give the correct view of the 
airplane to improve the lift-drag ratio of the vehicle at the expense of a system that produces 
thrust by the burning of fuel. 

The initial motivation for the use of active boundary layer control on aircraft 
surfaces was the growing need to conserve energy. The initial studies related to laminar 
flow control were led by NASA Langley and concentrated on the improvement in airframe 
L/D (see refs. [1] and [15]) although some attention was given to the arrangement of the 
propulsion system (see [16]). The range contribution of the propulsion system was not 
thoroughly analyzed, this will be the purpose of the ensuing discussion. In particular, it 
will be shown how the range is maximized by integrating the airframe boundary layer with 
the propulsion system in order to improve the airplane lift-drag ratio and the propulsive 
efficiency in combination. The term conventional aircraft will denote an aircraft powered 



by a conventional engine (air-breathing engine or propfan) and will be referenced as 
airplane A, while the non-conventional term applies to an aircraft powered by an air- 
breathing engine and a boundary layer thruster and it will be referenced as airplane C. 

The Brequet range equation assumes cruise flight at constant velocity V 0j with the 

thrust F equal to the drag D and the lift L equal to the weight W. This equation may be 
written as 

R = hii 0 L/Dln(^&) (5a) 

Wg 

where 

Tlo = *lp * Tit X Tin- (5b) 

This equation was originally derived assuming that the lifting and propulsion systems were 
separate systems. In the following analysis the fuel heating value, h , the thermal or cycle 
efficiency, T| t , and the transfer efficiency, % , are all assumed constant.* 

The propulsion efficiency , rip , the thrust per unit mass flow rate and the airplane 
L/D ratio will be discussed as they relate to a conventional aircraft. Finally, the product of 
the propulsive efficiency and the lift - drag ratio will be dealt with as it applies to an active 
boundary layer control airplane. The air-fuel ratio and the weight of the suction system- 
thruster configuration will be neglected in the present analysis. 

4.3 Propulsion Efficiency 

The propulsive efficiency, r) p , is defined as the useful work output per second of 
the propulsion system, FV 0 , divided by the energy input to the jet 1/2 m ( V e 2 - Vj 2 ). The 
latter can, in the limit, be the shaft input to the propeller or the fan of a high bypass ratio 
fan. Figure 9 is a schematic of the power flow for an aircraft employing a propfan and a 
boundary layer thruster. 

The propulsive efficiency is given by, 


* note that T| t is constant at the value it would have if the inlet for the power producing unit (not 
the thrust producer ) were in the freestream (see Fig. 10 of Ref. 15). 



(6) 


T, P = FVo - m(V e -'Vi)Vo 
E ( jd ) 1/2 rh ( V| - V? ) 

2 

v e / Vo + Vi / Vo 


For a propeller in the freestream Vj / V 0 = 1.0 and in the limit of V e / V 0 = 1.0 eq. (6) 
reveals that for these conditions ri p = 1 .0. In a boundary layer thruster using the airframe 
boundary layer the intake drag is charged to the airframe by virtue of the frictional process 
occuring in the boundary layer. In the limit of zero inlet momentum and V e / V 0 = 1 .0 , it 
is found that r| p = 2.0. This value of the propulsive efficiency is achieved considering an 
ideal intake system with zero losses. The attainment of zero velocity of boundary layer air 
is possible at the expense of friction. In the case of a flat plate this may be realized in the 
limit of an infinitely long plate as has been discussed by Pfenninger [17], In such a case 
the propulsive efficiency may exceed unity. It is also worth noting that if 


'Ipi 


_ if) Vj Vq _ 2 

l/2mV? V i/ V 0 


and 

'Hpe 


m V e Vq 
1/2 rh V| 


2 

Ve/Vo 


(7) 

( 8 ) 


then from eq. (6) 


Up — * 

2 / T^pe + 2 / t| p j 


(9) 


so that for t| p = 1.0; Upe = r| pi = 2.0, for example. 

This is the result described on page 16 in appendix B of ref. 16, although not 
directly in these terms. Equation (9) and the last form of eq. (6) imply that the exit and inlet 
terms add together, the next to the last form of eq. (6) shows that the inlet term, Vj , 
subtracts from the thrust and energy input to the engine. 

V. 

4.4 Thrust per mass flow rate of air, F / m , through the propulsor 

One of the distinct advantages of suction from the boundary layer in relation to the 
freestream is that it provides a larger thrust per pound of ingested flow. Considering the 
basic thrust equation 
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( 10 ) 


F = Ve . Vj_ 
m Vo Vo Vo 

for which Vj / V 0 = 1 , in the limit , for a propeller or a high by-pass ratio fan ( propfan ), it 
can be seen from Table 1 . that for the same value of V e / V 0 and in the limit for a boundary 
layer control system V) / V 0 = 0, the thrust per unit mass of ingested flow for the boundary 
layer thruster can be many times that of a conventional propeller system or an air-breathing 
engine. 


4.5 Airframe Maximum Lift-Drag Ratio 

The airplane lift-drag ratio in eq.(5) assumes that the propulsion system and the 
airframe are separate. This case will be considered initially. For this case 


L = 

D C Do + K C l 


(ID 


where the coefficients are based on the wing plan form area; i.e., C L = L/qoS where q 0 = 
1/2 pV 0 2 and K * 1 / it(AR)e where AR is the wing aspect ratio and e is the lifting 
efficiency. The term KCl 2 is the induced drag coefficient, Cp;. The L/D ratio may be 
maximized for constant K and Cqq as follows 


d(L/D) _ Q 

(12) 

<1Cl 

Here 


C Di s K C L op t = C °0 

(13) 

where 


CL ° P t - V K 

(14) 

It is easy to show that 


(L/D) max = CLopt = J- /, / — 1 — 

(15) 

2 C Do 2 v K Cd 0 


i.e., the maximum lift-drag ratio occurs where the induced drag C Di equals the drag at zero 
lift or the parasite drag, . 


4.6 Maximizing the Range by integrating the propulsion system with the airframe 

The conventional airplane or reference airplane A has a common propulsion system 
which overcomes the parasite drag and the induced drag and is completely separate from 
the airframe. The non-conventional airplane, airplane C, has separate thrusters for the 
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boundary layer and induced drag systems. In current practice the range is maximized by 
maximizing the L/D ratio and the propulsive efficiency independently . It will be shown 
however, that the use of an integrated system results in a smaller value of the lift coefficient 
which implies that the airplane altitude ^considerably reduced. The inherent advantage in 
flying at a higher altitude is an increase in the thermal efficiency of the air-breathing engine; 
however, flying at lower altitude results in an increase in the inlet flow rate to both the 
boundary layer thruster and the air-breathing engine for the same corrected flow rate and 
the same 'match' point (see [17]). 

The net thrust produced by airplane C is given by 

F c = rri c (V e - V 0 ) + m' bl V c (16) 


and the power P c of this aircraft is 

P c = l/2ni c ( V| - vj) + l/2m‘ bl V| (17) 


In the subsequent analysis it is assumed that the thrust provided by the main engine 
overcomes the induced drag while that provided by the boundary layer thruster overcomes 
the parasitic drag. The propulsive efficiency for this system may be written as 


- useful energy out _ D Vp 
P jet energy in Dp Vp + Dj Vp 
^pbl hpi 

_ Cp 0 + K C b 
Cop + Kcj 

hpbl hpi 

From the above equation and the range equation, eq. (5), 

I Con + K C? r. 

R = const. T| p — = 2 k LL. 


(18) 


D 


Cop + KC|: C D p + K C l 

hpbl hpi 

c L 


■ const. 


C Dp , Cp, 

"Hpbl hpi 


• const. 
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The optimum lift coefficient, the value that maximizes the product of ti L/D is found from 


d (t| p L/D) _ ^ 
dC L 
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( 20 ) 


where the resulting lift coefficient is 


c, = yi/£20 Uhl 
L °P l V K *lbl 

Thus the ratio of the range of the non-conventional airplane to the conventional airplane is 



which for (Vj/Vq)^ = 1-0, (V e /Vo)i = 1.25, Ti pb i = 2.0 and Ti pi = 0.89 the range ratio is 
1.5. 


From the preceding analysis and subject to the underlying assumptions it is seen 
that the range of the non-conventional aircraft employing a boundary layer thruster is 
approximately 50% greater than that of the conventional. It is estimated that for a parasite 
drag coefficient of 0.015 the change in altitude is about 9,000 ft. for a conventional aircraft 
flying at 40,000 ft. with a wing aspect ratio of 7.0 and a lifting efficiency factor of 0.92. 

The preceding analysis was performed without recourse to the flow pattern over the 
aircraft surfaces; the presumption being made that the boundary layer was turbulent and 
that it remained unchanged for the case of the non-conventional aircraft. The main purpose 
of bleeding the boundary layer off the airframe is to convert it to a laminar layer. This 
conversion results in a considerable reduction in the friction coefficient as may be seen 
from figs. 10 and 1 1 for the flow over a flat plate (see [4]). At a Reynolds number of Re = 
10 7 the boundary layer is usually turbulent and the ratio of the turbulent skin friction 
coefficient to the laminar one is 8 : 1. If suction is applied such that the friction coefficient 
with boundary layer suction is twice the laminar value with no suction, the ratio of 
maximum L/D for the bleed case to the non-bleed case is 2.0 from eq. (15). The ratio of 
the corresponding lift coefficients becomes 1/2 which would indicate an optimum cruise 
altitude of 16,000 ft. lower than for the conventional aircraft with a turbulent boundary 
layer for the same wing loading W/S. 


5. Conclusion,- 

A brief study on the use of suction air and the associated propulsive benefits has 
been made. It has been shown that the intake mass flow from the boundary layer offers 
more thrust per pound of ingested air than the freestream. In particular, it was shown that 
in the re-acceleration of boundary layer air to freestream velocity the concept of a 
propulsive efficiency greater than one is not far-fetched. It has also been demonstrated that 
the propulsion of boundary layer air via an auxiliary propulsor entails a different 
perspective of the Brequet range equation; specifically, one that deals with an integrated 
system and which incorporates the propulsive efficiency of the thruster and the lift-drag 
ratio in determining the aircraft range at subsonic speeds. With this view in mind, the 
possibility of flying at lower altitudes to attain the desired range is realized. 
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A-l Derivation of the Payload equation 

The payload equation of a propulsive system may be derived as follows , this is 
done for a conventional aircraft, i.e. an aircraft not employing a boundary layer thruster, 
for an airplane at cruise. 

The range equation may be written in alternative form as 
x = 863.5 y^logio(^j) (miles) (A.l) 

The initial overall weight of the aircraft at takeoff, Wj, is composed of the weight of 

fuel, w f , the structural weight, w sl , and payload weight, w L , 

> 

Wi = Wf + w s[ + w L (A. 2) 
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Let us initially consider that the takeoff weight is equivalent to the cruise weight at a 
certain altitude at zero range ( this is not entirely correct and varies with the type of aircraft 
and its takeoff speed.)* 

Substituting the expression for the payload into the range equation. 


_cx_ 


863.5 Ti^ 


= logJ 


1 | w st , W L - 
Wf Wf 
W st , W L 
Wf Wf 


(A. 3) 
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wl _ 

Wf 
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863.5 n 2- 


Cd 


W 5t 

Wf 


The payload to overall initial mass ratio is given as 


wl Wf 

W 1 1 + Wf/Wf + W st /Wf 


Now let 

[c/863.5 r\^]= p 

then 


w L _ (l + W st /Wf) Wst .. .. 

Wf o Wf -h (A. 5) 

10 * 


* Note that if we were considering a non-conventional aircraft the weight of 
the suction system-thruster configuration must be included as part of the 
structural mass of the aircraft. The lubricating oil is included in the total fuel 
consumption. 



The ratio w^/wj = 1 represents the total weight at zero range ( or take-off weight). 
The diagram shown below is a typical payload-range diagram of aircraft capability. The 
diagonal range line is limited either by maximum takeoff weight or by the fuel capacity of 
the aircraftThe outer curve represents the projected curve for the non-conventional aircraft 
employing laminar flow control with propulsion and the possibility of maximizing the 
payload capability over any given range of the aircraft. A lower value of the parameter p 
indicates an increase in the payload as suggested by equation (3). 


Payload 


Maximum Take-off 



A-2.1 The propulsive efficiency 


rj in the above expression is the propulsive efficiency of the entire propulsive 
device. If we have an aircraft powered by a single engine then i\ is the propulsive 
efficiency of that engine. Now suppose we have the same aircraft equipped with a 
boundary layer thruster ( a nonconventional engine) and an air-breathing engine( a 
conventional engine) what would be the net propulsive efficiency? 


First, let us consider an aircraft with two conventional engines, the propulsive 
efficiency is given by ^ 


t 1p = 


Pj total thrust power 

P e total power output of engine(s) 


P t = (Fi+F 2 )U , P e = P T +(K.E)e 
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Let us suppose that the two engines have different jet velocities. 


(K.E)c =(K.E) cl +{K.E) c2 

where 

(K.E)ei =^(l+fl )v2j = ^(l+fiJ^jj-U) 2 

(K.E)e 2 = ^(l+f 2 )v 2 2 = Uhl (1 +f 2 )( Vj2 -U)2 

or we can obtain rip directly from the ratio of useful output over input and neglecting the 
fuel-air ratios 

_ tti a i (Vji - U) u + Iti a2 ( Vj2 - U) U 
P dial ( v 2 , - u 2 ) /2 + rii a2 (v 2 2 - U 2 )/2 


The above analysis, however, is done with undue regard to the airframe 
aerodynamics. In other words, what effect do two jets emanating at different velocities 
have on the drag. Another question in connection with the airframe is how would the 
engines be mounted such that we have two different mass flow rates issuing into the 
engines. 

One question in connection with the term efficiency is this: In dealing with the 
boundary layer thruster can we speak of both a thermal efficiency and a propulsive 
efficiency of the boundary layer thruster? We may answer this question in the following 
manner 

If we are looking at an air-breathing engine we may define both a thermal and a 
propulsive efficiency. The thermal efficiency as a result of the diffuser or inlet - 
compressor-burner or combustion chamber - turbine - configuration, and the propulsive 
efficiency if a propulsive device is attached to this main section. The boundary layer 
thruster on the other hand is a propulsive device'and can only possess a propulsive 
efficiency. 

The thermal efficiency of a device is defined by the work done over the net quantity 
of heat added. This parameter cannot apply to a boundary layer thruster since there is no 
net heat added ( the entire process may be considered adiabatic) and the work may or may 
not exist depending on the particular thrusting device employed. If the device is an 
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accelerating device such as a nozzle the net work done is zero. On the other hand for a 
suction compressor the net work done is not zero. 

A-2.2 On the total drag or total thrust 

When dealing with the total drag of an aircraft why are we concerned with drag 
values of the components of the entire aircraft? One major contribution of boundary layer 
suction is the considerable reduction of drag over that of a conventional aircraft; as such, an 
accurate prediction of drag is important for estimating the benefits of boundary layer 
suction. Let us first consider the estimation of drag of a conventional aircraft fitted with a 
single engine. 

(i) method A ( direct approach): 

Suppose we have an aircraft fitted with a single engine traveling at certain altitude 
with a given speed. With these values we can determine by a thermodynamic analysis of 
the engine cycle the thrust generated. A question at this point is how are we sure that the 
thrust generated is the correct thrust needed to propel the airplane at the given speed? If we 
can determine the necessary thrust then we know that this is the total reaction of the air 
upon the aircraft. 

(ii) method B ( indirect approach): 

This approach relies on determining the drag directly from knowledge of the drag 
values of the components of the aircraft 

If we consider the use of a boundary layer thruster in conjunction with the air- 
breathing engine, then for suction over the wing the entire wing drag (neglecting 
compressibility effects) is the sum of wake drag, suction drag and induced drag. This must 
be less than the profile drag of an unsucked wing of the same planform to justify the use of 
suction. The maximum lift varies considerably with suction since for suction over the top 
surface of the wing the flow is accelerated with the consequent reduction of pressure over 
that surface. * 
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Appendix B 

B-l. Thermodynamic relations used in modeling the operation of the boundary layer 
thrusting device. 

The following section deals with the governing equations used in modeling the 
boundary layer thruster based on the simplified schematics shown below. The process of 
the deceleration of air from freestream conditions to the boundary layer edge is considered 
to be an isentropic process. The collection surface may be taken as a flat plate with equality 
of the static pressure at the surface with the local pressure at the boundary layer edge, 
process 1-2 from the collection surface to the entrance of the suction compressor inlet is 
considered to be an adiabatic process with input of approporiate pressure drop values into 
the computer program. The compression process is also taken as adiabatic. The relevant 
equations appropriate to the individual processes are given below. 


FreeS ir earn 
T~ P- M~ 


M- 


Boundary Layer Edge 

\ 
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LFC Surface 



Figure B.l - Schematic of Boundary Layer Thrusting Device 


a. Isentropic flow of air from freestream condtions to local values at the boundary 
layer edge. 

Certain conventional defintions are given below : 
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M -=ar= 


VyRT: 


hxoo - hoo + 


(B.l) 


K 

Cp T Xoo — Cp T„, -i — y RT„ 


1^=I + ^M 2 
T„ 2 “ 

for an isentropic process Pv Y = constant, therefore 

JV. ITtJ^ 

P„ \ j 

by definition the dimensionless pressure coefficient is given by 


(B.2) 


(B.3) 


C P = 


P-Po. 


1/2 P-X 

the local pressure is Pj and 

^ _ P!-P~ 

Cpi = 7 

i n p. vi 

from which follows that 


(B.4) 


The magnitude of the adiabatic wall temperature T aw relative to the one-dimensional mean 
values of the static temperature T and the stagnation temperature T 0 is expressed by the 
recovery factor R f , where 

R f = Ir 

T 0 -T - 

The average stagnation temperature in the external flow is taken as the freestream total 
temperature. The assumption is made that the total pressure of the flow at station 1 ( the 
LFC surface ) is the same as the external static pf^ssure i.e.; P X1 = Pj . The local 
temperature and pressure values at the boundary layer edge are related by the following 
equation 

T ,^= T -(^jr= T .(^)T 

b. adiabatic process in suction duct : process 1-2 
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The adiabatic wall temperature is taken as the stagnation temperature at the LFC surface T aw 
= T t1 and since process 1-2 is adiabatic, T aw - Tti - t T 2- Using the definition of the 
recovery factor given above, the expression for the total temperature at the entrance to the 
compressor T T2 may be found as follows. 

(Tpoo - Tjocal) Rf = Tt 2 " Ti oca i 
Tj2 = T] oca l + (T Too - T] 0 cal) VPF 

where for a laminar boundary layer the Prandtl number is approximately 0.7. 


Tt2 = 




= T (Si.fr 

00 \Poo| Y 

= T -(# +VPF ^--T.vpF(£f)T 

= T_(l-fl^a)T+T_VFr^ 

\ Fo °/ loo 

J-l 

( 1 - VPr I^-jy + VPr 


(B.5) 


The relative pressure loss between the collection surface and the inlet to the compressor is 


PTl - PT2 _ Pi - PT2 
Pti Pi 

PT2 _ Pi /, (pi -PT2) ) | 

Poo Pool 1 Pi ) 

c. adiabatic compression process across suction compressor : process 2-3 A 

This process may be represented on the h,T- s diagram with the isentropic path shown as 

2-3A' . 
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Figure B.2 - Compression Path of suction flow air across adiabatic compressor. 
The adiabatic compressor efficiency is defined as, 


_ hr3A' - h T 2 _ T ( T3 a 7T T2 - 1 ) 
lc hT3A-h-n T2 lTT3A-T-ni 
rearranging, 


Tt 3A = T-n 


J-(T 3 a/Tt2-1) + 

tic 


1 


but 

Tt 3A' _ ( PT3A' V^r _ / PT3A V~~ 
Tt2 V PT2 I 1 \ PT2 / ' 


therefore, 
Tt3A = Tt2 



+ 


1 


The steady flow energy equation per unit mass is given by. 


(B.7) 


(B.8) 


V 2 


h2 + “^ + gZ2 ± qext + W c = h 3A + + gZ 3A 


(B.9) 


where q ext is the external heat loss or gain per unit mass of air flow into the suction 
compressor, and w c is the work added to the system per unit mass flow of air. Neglecting 
potential energy changes and since the process is adiabatic the above equation becomes 
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y2 y 2 

h2 + -*■ + w c = h 3A + -&■ 

hT2 + w c = hl3A 

substituting for Tj 3 A , w c is given by * 



d. process 3 A - 3 : The air flow form the exit of the suction compressor to the entrance 
to the C-D nozzle is considered to be adiabatic. The total pressure loss Ptca ' Pt3 / Pt3A * s 
taken as an input into the code. 

e. process 3-e : The nozzle expansion process may be represented on the h-s diagram 
as shown in figure 3.C. 




Tt3 = Tie 




Figure B.3 - Expansion Process Through C-D Nozzle 
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The adiabatic nozzle efficiency is defined as 


_ h-n - h e _ Cp (Tt 3 - T e ) 
hj3 - h e ' T-n - T e - 


where. 


Tt3-T c = 


yL 

2 Cp 


therefore, 

_ v| / 2 Cp _ Ve / 2 Cp 
^ Tt 3 - T e ' (1-T e */TTJ)TTJ 


(B.ll) 


but, 

lel = (P£l)Y 
T-n VPT3/ ' 

here P Te’ = P e = P oo and the expression for the exit velocity becomes, 

V.=V 2c p’1" T T3('-(^) . 

the ratio of the nozzle thrust power to the compression power is given by V e V M / w c . The 
Mach number in terms of the exit temperature is given by M e = V e / ( yR T e ) 1/2 . The exit 
nozzle area may be found as follows. 


Hi 

A 


= pV = p V/RT = 


P ,/Y,/T7 1 
VyRT V R V T VTq 


= VI^ M V i+I f M2 <B - 13) 

p may be eliminated using 
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M 


= Jl- Po- 

V R fTo 


in terms of the exit area and properties at the exit A e is 


Jy + 0 

1 +^r-M 2 2(Y- 0 


Ae = - 


lhair VRT-n 


-( 7 + 1 ) 


PreVYMe 1 + 


where P Te is the total pressure at the exit of the nozzle given by. 


(B.14) 


■»- p - 

The numerical code implementing the above equations is given Appendix C together with 
the numerical values obtained for the input parameters given below. Figures B.4 and B.5 
are the performance plots of the boundary layer thrusting device based on these numerical 
values. The ideal curve, curve A of Figure B.4, represents a zero total pressure drop in the 
duct leading from the LFC surface to the inlet of the suction compressor and for that 
leading from the exit plane of the compressor to the inlet plane of the C-D nozzle with ideal 
efficiencies of 1 for the compressor and the nozzle. Curve B represents a zero drop in total 
pressure with efficiencies of 0.8 and 0.98 for the compressor and nozzle respectively. It is 
seen that the ideal curve models very closely eq. (6) with Vj = 0. Figure B.5 shows that 
the optimum value of the exit velocity is the freestream value for the lowest coefficient of 
incremental drag as derived in the following section. 
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Fig. B.4 - Thermodynamic Performance of LFC Pumping System 
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Fig. 6. 5 - variation of incremental drag per suction flow rate per suction area ratio 


The input flow parameters representing the non-ideal cases are as follows: 




M eo = 2.2 

T„ = 216.65 K at 60,000 ft. 
p„ = 0.1 1532 kg/m 3 
c p i = 0.0 

Cp = 1011.5 kJ/ kg. k 

R f = 0.7 
Tin = 0.98 
Tic = 0.8 

R = 287 kJ/kg. k 
y= 1.4 

PSFC = 0.3308 lbm / hr - hp 
TSFC= 1.28 lbm/hr-lbf 

PTC*' PTC = 0.05 


Pl 'PTC = 0.00 , 0.10 , 0.20 , 0.30 


C.2 - On the minimization of overall fuel flow. 

The total amount of fuel required for the operation of an air vehicle coupled with a single ( 
or multiple) conventional air-breathing engine and a boundary layer thruster is the sum of 
of the fuel flow required by the air-breathing engine plus the additional fuel required to 
provide auxiliary shaft power to the LFC compressor, 

(tflf)lotaI = (tilf)main engine + (^f^LFC pump 

By definition, the thrust - specific - fuel consumption TSFC is 
TSFC — en g* n e 


where x main engine is the thrust provided by the main engine 
^main engine” total drag - T n ozz]e — Diota] -T nozz ] e 

where T nozz i e is the thrust provided by the LFC thruster. The power-specific-fuel - 
consumption is given by : 


PSFC pump 

Ppump 

The pump efficiency ri ps is given by 

Tnozzie 

Tips = -5 

r pump 
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• » 


_ TlpsP pump 

^nozzle ” y 

* oo 

assuming a transmission efficiency of r|tr= 1.0 in the conversion of shaft power from the 
turbine to the auxiliary compressor , the total fuel flow rate required is, 

(dlf)total = TSFC (Dtota] - Xnozzle) + PSFC Ppump 

= TSFC (D tota] - T nozz]c ) + Ppump = TSFC(D e ff.) 

where the effective drag is given by 


D eff . = D lotal - T| ps -Ppi + P pump 

By definition the lift coefficient is given by C L = 2 L / p^V^S = L / qS. The effective drag 
may then be written as 


D rff . _ TSFC V. 


pump 


CDeff. „ c „ c 

(j S (| s 

The total drag coefficient is C Dlolal = D tolal / qS therefore. 


CDeff. = ^Dtotal + ( 


I TSFC V„ 

_ Cl q S _ 


pump 


qS 


Cl. 


Deff. CDeff. q S 


CDtotal + (ygpQ V oo - Xl p s| 


r pump 

SqV„ 


the term 


( PSFC v „ \ P pump 
'TSFC 00 lps / S q V„ 

is defined as the fuel-equivalent drag increment AC d . The pump power in the above 
expression P pump is the suction mass flow rate multiplied by the suction compressor work 
done per unit mass flow rate as provided by eq.(10) above. The results presented for the 
parameter AC d S / (pu) s Sj are shown in figure B.5, which is the incremental drag as 
defined above per suction flow rate per suction area ratio, which takes into account the 
equivalent suction area S] and the volume flow rate of ingested boundary layer air. 
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Figure 1 . - Drag coefficients for the flat plate at zero incidence 
with uniform suction. 

Curves (1), (2) and (3) refer to tk> suction 

(1) Uminar 

(2) transition from Uminar to turbulent 

(3) fully turbulent 
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Figure 2. - Comparison of ouising range. 
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Figure 5. - Effect of HLFC on Wing Profile Drag 
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Fig. 8 - Schematic of re-acceleration of ingested flow at the 
trailing edge of a flat plate to frccstrcam velocity 
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Fig. 9. - Power Flow Diagram for a Jet Airplane Employing a Boundary Layer Thruster. 


3 9 





4*1 
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Figure 10. - Coefficient of skin friction of a flat plate at zero 
incidence. Optimum suction denotes smallest volume coefficient 
required to maintain laminar flow. 



Figure 11 . - Relative saving in drag on flat plate at zero 
incidence with suction maintaining laminar flow at optimum 
suction 
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Computer Code and Numerical Results 



TUB PRC-uRA;* SOLVES FOR T-i I ABLE ° A R A!i BTE » S OF THE THcSMO- 
DYNm.IC RELATIONS APPLIED TO THE LFC ^UMFInG SYSTEM. 

VARIABLE DIRECTORY : 

MINF : THE FREE- STREAM MACH a 

TTTF : TrTE FR77-ETR : TTX? FRATUR T 

P I N F : THE FREE-STRsA", e R ISSUES 

CcL^C : THE RELATIVE i iO? IN TOTAL PRESSURE BETWEEN STATIONS 
1 t 2 

D”EL”P3 : — Pm EX ATT VE UTC? IN TOTAL ~PF?FS'SURT _ ?rrW rE N "ST ATI C N S " 

3 A Li 

CPI : PRESSURE COEFFICIENT 

RP3 : THE RATIO OF TOTAL PRESSURE AT STATION 3 TO THAT AT 

FRFr^TTRTT^ 

ETAC : T^E COMPRESSOR EFFICIENCY 
E T A N : T H - NOZZLE EFFICIENCY 
G A v , A : KaTIC CF SPECIFIC h £ •• T S 

Tz : TTOZnrC - TXTT*TETOrTTT 

TT3 : TOTAL TEMPERATURE AT STATION 3 
PT3 : TOTAL PRESSURE AT STATION 3 
FNOZ • N 0 1 Z u E TriRUiT 

FR FThF-FT JT^IZ NyX.SR 

MAI* : F 1 3 j F L 0 a RATE L r A I F. 

A E : NOZZLE EXIT AREA 

PCOHP : COMPRESSOR *'0e< p = p ft IT MASS OF AIR 

Fcrar^~r"R att^tif^otzue ~ ttfjtt ? ove r hpt rtci sit mass - wars tc ~ 

PCOVP 

ME : EXIT MACH NUMBER 


REAL MI NF, R P3 ( 33)/ ME/M A IR/ POWER (33) /VELP (33) /DELP 2(4) 

REAL PI (<*/33)/PZ(4/33) 

DATA R°3/j.094 / 0.095 /. 093/ .Q9j /« 0 9 9/3. 10 / »1 c/. 14/G .1 5/Q. 2/ .2 5/ Q.3Q 
c* / • 3 1 / 0 • 4/ . 3 / • 5 5/*0/aDj/« 7 / D . 7 ) / j > o / j « 3 3/ 0 . v / 0 . 9 3 / 1 . j / 1 . 3 /7 • u / c . 3 / 
6.3. 0/3. 5/4.0 /A. 5/ 5.0/ 

DATA DELPZ/ J.GU/C.1 /j.^/0.3/ 

INPUT TH£ CONSTANTS OF TnE P-0:- LEX: 

£ R 1 T t TT/^T^T^'P J T THE CONSTANTS OF TW P»OSTe'fls T 

WRITEtS/^I'MINF/TlNF/ROINF/GAMA/R/PR/CPI/ETAC/ETAN' 

R£AD(5/*)MINF/TIsF/R0lNF,GAMA/R/Pft/CP1/ETAC/ETAN 

»JRIT£C7/*) , MINF/TINF/ROINF/3AMA/R/PR/CP1/ETAC/ETAN , 

• AITE( 7 7 R I'FTKTTTTK V / R C T uT7 1 IK7T7X7T r P. ,'Z /7TAT / E T An — 

PRINT*/* RP3 IS * 

PRINT*/RP3 
DELP3 = 0.3C 

”*~mTT5 

SAMAR = (GAMA -D/SAMA I 

MAI R * 1 .0 


ORKtfS^AL PAGt !S 
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FOR EQUATIONS 1-* ; 

TIFT = R 0 1 M r*’T*TT7rF 

VINF = MINr *3iRT<uA>!4*R*T INF ) 

VAR1 = 1 + <6A"A - 1 > /2*K INr** 2 
VAR2 = VAR 1 ♦ » ( 1 / 3 A w A R ) 

TTC ^r *“TC n TZ RTTAvi 

TT2 = T I N F * ( ( 1 - 5 3 R T ( P R ) )'* V h R 3 * * j A A \ 
«fttlT£(6s*) 'VlNf/VARI^VARZ/VAPi^TT 2' 
xKIT£(o^»)VINF^VAP . 1 /VAP. 2*VAR3*TT? 

DO 2 K = 1,1 

WRITc(7/*) , DELP2^DELP3^CP^oA'’AR/rfAIR' 
tiRITc(7s*)D£LP2(K)sDeLP3sCFsSAMARs*. AIR 
DO 1 I *1,N 


+ S:RT(?R)*VA-t1) 


s RP 3 C I ) *7 


H t i 


FOR cm J AT I ON S 7-13 


VAR3 = TA'R 3*(1-DcLPt(<)) 

V AR 6 = 1/VAR4 

V AR7 = VAR4/C1 -0ELP3) 
i/AS'-t = wa ST/VAR 5 


f n < \ Tf XI V u w r ^ j 

VAR ? ~ VAR7/VAF. 5 

TT = TT I*(1-CT-(w l A-;i)*'*oA' | A D7 "£Ta"C T 

PA;j 1 = 2* ETAS *C i3 *TT3*(1 - (VAF.£i**IAPA')) 
PRInT*/'PA?.1 = ' , ^ A 4 1 

_ v i = s:rt( 3 aad 


P C 0 M ? = CP* T 

P0«£R<I) = VE*VINF/PCG*P 
VElR(I) = VE/VINF 


rr a c * < < vtrt**oXka ^'5 


■IHfcB-HfcfHh 

TE = T T 3 - VE**2/(2*C?) 


ME = VE/3iRT(3A\A*R*T£) 

PTE - PINF»< TT3/T£)»«(1 /^A^Ak ) 

C - = Cl + C GAMA-1 JTI*TZ * * 2 ) * * C -TTS vA' _ F r T77‘CZ *TGA FA-1 ) ) ) 

Ai. = /'i AI R*3 JRT(R*TT3 )/ CPTE* 3 j RT ( GAM A) •’it *Di.’») 


OUTPUT RESULTS • 


WRIT £(*✓*) * VA 3 ‘»^V^R5^VArti^VA-<7/VAR:zTT3^' > Ar'1 * 
*UTE(:5,*>VA?4,V*4G*VAR'j,VAP7rV4R3/'TT3* 1> AR1 
*3.1 T £ C o_c *) ' V 3 f P £ 0 •: P / V E l 3 , P 0 « cR/T E^IEOE'wAE' 
77TT77 /»,, i ^ .JZ iVTosT Tijn / r rrrc-rr ~vr— 


— ' V-/ p Cu-:P/ VE l p / P 0 b E/ f 1E/DE'wAE* 

1 k t x pcoWp/vel « rry/p o * £ * ( i > ✓ t et^e/dtnttt 

CONTIN'Jc 

?3sI?4^ 4 y > 4 v - L ?5 ^^overcd^isi^m) 

FORMA T C 1 / 2 F 1 u . 5,4X,F10.5> 

4A7TtT«r71TJ 


^ Ts -S$?2j) (I/RP3(I>,yELR(I),P0«!eR(I>, 1 = 1 r O 

2?-^IS! !'£ x ' 1 *,3X,* kf-O',?*,' VELP.'x 11 X,» POWER 

0 K A T ( ' , r2XxT4,3X.F7. 4 < r5X^F1J. 5»5 X J >F1C.5) 


COKTIKJS 

= < 7' ? C> (aP3<I ) ,°1C1,I),P2<1,i),Pi (2,1), P2( 2, I),P1 (3,1), 
5P2(a#I)/P1(4zI)^P2(4^I)^I=1^N) 


to 


V 


FORM AT < * '^Fi.4/2XrFG.4, 
SFc."»^2X/Fo.4r2X / F : . 4) 

STOP 

END 


;X^Ft.4^2X^Fi.4^2XxFw.4^2X/Fo.4^2X/ 
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Page #1 


"lfcpuir\p.data2" 


Tuesday, June 12 3:40 PM 1990 


Pt 3 / Ptoo 

0.094000 
0.095000 
0.096000 
0.098000 
0.099000 
0.10000 
0.12000 
0.14000 
0.15000 
0.20000 
0.25000 
0 . 30000 
0.35000 
0 .40000 
0.50000 
0.55000 
0.60000 
0.65000 
0.70000 
0.75000 
0.80000 
0.85000 
0 . 90000 
0.95000 
1.0000 

1.5000 

2.0000 

2.5000 

3.0000 

3.5000 

4 . 0000 

4.5000 

5.0000 


(Ve / V«) 1 (Power ratio) 1 


0.052300 

2.7391 

0.091900 

4.0429 

0.11880 

4.5133 

0.15930 

* 4.7630 

0 . 17590 

4.7599 

0.19100 

4.7225 

0.37560 

3.5592 

0.48550 

2.9738 

0.52910 

2.7858 

0 . 69110 

2.2511 

0.80380 

1.9863 

0.89120 

1 . 8215 

0 . 96280 

1,7063 

1.0237 

1.6198 

1.1236 

1.4960 

1.1658 

1.4495 

1.2042 

1.4098 

1.2393 

1.3753 

1.2718 

1.3450 

1.3020 

1.3181 

1.3303 

1.2939 

1.3568 

1.2720 

1.3817 

1.2521 

1 .4054 

1.2338 

1 . 4278 

1 .2170 

1.6058 

1.0985 

1.7336 

1 . 0270 

1.8339 

0. 97730 

1.9167 

0 . 93970 

1.9875 

0.90980 

2 . 0495 

0.88520 

2.1046 

9 0.86440 

2.1544 

0.84640 

0 . 0000 

0.0000 


(Ve /V°o)2 

(Power ratio) 2 

0.053300 

0.95830 

0.093600 

1.5767 

0 . 12110 

1.9188 

0.16220 

2.3041 

0.17920 

2.4200 

0.19460 

2.5060 

0.38260 

2 . 6436 

0.49440 

2.4212 

0.53880 

2.3254 

0.70360 

1.9967 

0.81820 

1.8060 

0.90690 

1.6790 

0 . 97970 

1.5866 

1.0415 

1 . 5154 

1.1430 

1.4112 

1.1858 

1.3712 

1.2248 

1.3368 

1.2605 

1.3068 

1.2935 

1.2801 

1.3242 

1.2564 

1.3528 

1.2349 

1.3797 

1.2154 

1.4051 

1.1976 

1.4291 

1.1812 

1.4518 

1 .1661 

1.6326 

1.0584 

1.7622 

0 . 99250 

1.8640 

0 . 94630 

1.9481 

0.91110 

2.0199 

0.88310 

2.0828 

0.85990 

2.1387 

0.84030 

2.1892 

0.82330 

0.0000 

0.0000 
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"If cpump . data2" 


Tuesday, June 12 3:45 PM 1990 


(Ve / V°o) 3 

(Power ratio) 3 

(Ve / Voo) 4 

(Power ratio) 4 

0.054500 

0.55660 

0 .055700 

0.37810 

0.095600 

0.94010 

0.097900 

0.64630 

0.12360 

1.1712 

^0 . 12650 

0.81410 

0.16560 

*1 . 4645 

*0.16950 

1.0384 

0.18290 

1.5655 

0.18720 

1.1200 

0.19860 

1.6474 

0.20330 

1.1887 

0.39040 

2.0556 

0.39950 

1.6435 

0.50450 

2.0063 

0.51610 

1.6811 

0.54980 

1.9641 

0.56240 

1.6710 

0.71770 

1.7733 

0.73390 

1.5742 

0.83440 

1.6399 

0.85300 

1.4852 

0.92470 

1.5440 

0.94520 

1.4151 

0.99880 

1.4712 

1.0208 

1.3592 

1.0617 

1.4136 

1.0849 

1.3135 

1.1649 

1.3269 

1 .1903 

1.2426 

1.2086 

1.2930 

1.2347 

1.2143 

1.2482 

1.2636 

1.2751 

1.1895 

1.2845 

1.2376 

1.3122 

1.1674 

1.3181 

1.2145 

1.3464 

1.1476 

1.3493 

1.1937 

1.3782 

1.1296 

1.3784 

1.1749 

1.4079 

1.1133 

1 . 4058 

1.1577 

1.4358 

1.0983 

1.4316 

1.1419 

1 . 4621 

1,0845 

1.4560 

1.1273 

1.4869 

1.0717 

1.4791 

1.1138 

1.5105 

1.0598 

1.6629 

1.0167 

1.6978 

0.97300 

1.7947 

0 . 95640 

1.8321 

0.91820 

1.8981 

0 . 91370 

1.9375 

0.87910 

1.9836 

0.88100 

2.0246 

0.84900 

2.0567 

0.85490 

2.0991 

0.82470 

2.1206 

0.83320 

2.1642 

0.80450 

2.1775 

0.81470 

2.2221 

0.78730 

2.2288 

0.79880 

2.2744 

0.77240 

0.0000 

0 . 000 ® 

0.0000 

0.0000 




>age ll - "pump data 4" 


Tuesday, June 12 3:36 PM 1990 


(Ve / V°o) (Power ratio) 


0.052390 
0.091950 
0.11883 
0.15915 
0.17570 
0 .19072 
0.37282 
0.47961 
0.52162 
0 . 67579 
0.78147 
0 . 86257 
0.92860 
0 . 98441 
1.0755 
1.1139 
1.1486 
1 . 1804 
1.2098 
1.2370 
1.2624 
1.2862 
1.3087 
1.3299 
1.3500 
1.5092 
1 . 6230 
1.7121 
1.7857 
1.8485 
1 . 9035 
1.9524 


38 .167 
21.750 
■ 16.830 
* 12.566 
11.383 
10.487 
5.3645 
4.1700 
3.8342 
2.9595 
2.5593 
2.3187 
2.1538 
2.0317 
1.8595 
1.7955 
1.7412 
1.6943 
1.6532 
1.6168 
1.5843 
1.5549 
1.5283 
1.5039 
1.4815 
1.3252 
1.2323 
1.1681 
1.1200 
1.0819 
1.0507 
1.0244 


Appendix D 

Definitions and Aerodynamic Terminology. 

For a conventional airfoil the total resistance or drag is the sum of the profile drag and the 
induced drag. 

For an airfoil with a suction surface the total resistance or drag is the sum of the wake drag, 

the induced drag and the suction drag. The wake drag and the suction drag being the profile' 
drag of a suction surface. 

Profile drag : is due to the friction of the air along the sides of the airfoil 
induced drag : finds its origin in the circumstances that the appearance of the lift is 
accompanied by the creation of a definite flow pattern in the neighborhood and in the wake 
of the airfoil, which demands a continuous supply of energy. 

Suction drag, is not an actual physical drag acting to oppose the motion of the airfoil 
through the air.but a drag computed from suction power requirements; however, the 
suction drag may be considered as an actual physical drag when the suction flow rate is 
high enough that the effect of suction is felt by the external potential flow. This 
phenomena is then termed ’ the sink effect of suction '. 

Parasitic drag. In the complete structure of an aircraft are found various parts either of the 
structure or of the equipment which, like the body, take no pan in the development of lift 
and the drag of which may be grouped under the general head ’ parasitic '. This parasitic 
drag is made up of the fuselage, landing gear, tail surfaces, etc and of their interference with 
the wings and between themselves. Interference drag may be as high as the sum of the 
component of the drag of the component parts tested separately. 

Wake drag : The drag obtained by considering a control volume in the wake of the airfoil. 

This wake is an indication of a momentum deficit due to the presence of the body in a field 

of flow. A control surface taken upstream and downstream of the body is used in 

evaluating the profile drag from knowledge of the pressure and velocity distributions in the 
wake. 

Ram drag ; represents the loss of momentum associated with a conventional engine. 

For an aircraft with no suction surfaces: '\ 

total drag = profile drag + induced drag + parasite drag 
For an aircraft with suction surfaces: 

S ( profile drag of non-suction surfaces ) + induced drag + parasite drag + 
suction drag + ( wake drag of suction surfaces) 

Range. The distance that can be flown with a given amount of fuel. 


46 



Endurance: The hW . 

may be m0St] y a PP^cable to filter drpTaneTrail ^T' ^ Pm *^ of ^ormance 

planes ra.herfhancomrcercial subsonic transport 

ftopukive efficiency, = ratio of 10^"^^, i ^ mCn,uni of 8 ases + s™ of pressure forces 
TI.rase specific fuel consumption lsfc _ ^ '° ' 0,al e Power ouiput 

:-:::r=r=rz. 



